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ABSTRACT 

The  method  of  least  squares  is  applied  to  the  problem  of 
analyzing  a  decay  record  to  determine  the  damping  constants  for 
a  dry  friction- viscous  damped,  single-degree-of-freedom  system. 
Solution  of  the  set  of  non-linear  equations  which  yield  the  constants 
is  obtained  by  applying  the  Newton- Raphson  method  of  iteration. 

Sample  calculations  show  that  the  method  is  not  well- suited  for 
manual  computation. 

A  program  is  presented  for  calculating  the  damping  con¬ 
stants  by  means  of  a  digital  computer. 
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INTRODUCTION 


The  objective  of  this  work  was  to  develop  an  analytical  method 
for  analyzing  an  experimentally  obtained  decay  record  to  determine 
the  damping  constants  for  a  dry  friction-viscous  damped  s ingle- degree- 
of-freedom  system. 

A  semi- graph! cal,  trial  and  error  method  of  analysis  is  available;^ 
however,  as  pointed  out  by  Jacobsen  and  Ayre,*  a  precise  determination 
of  the  damping  constants  by  that  method  is  time  consuming.  Moreover, 
engineering  judgment  must  be  used  to  determine  when  the  constants  have 
been  calculated  with  sufficient  accuracy;  and,  the  semi- graphical  pro¬ 
cedure  may  lead  to  inconsistent  results. 

In  this  report  the  method  of  least  squares  is  applied  to  the 
problem  of  analyzing  a  vibration  trace.  The  method  yields  consistent 
results,  but  is  not  well- suited  for  manual  computation.  A  program  is 
presented  for  carrying  out  the  calculations  by  means  of  a  digital 
computer. 


ANALYTICAL  DEVELOPMENT 

The  present  study  is  limited  to  linear  systems  in  which  the 
damping  forces  are  due  to  a  combination  of  viscous  friction  and  dry 
friction.  The  objective  of  the  study  is  to  determine  the  damping 
constants  from  an  experimentally  obtained  decay  curve.  As  Indicated 
in  Figure  1,  the  experimental  curve  for  such  a  system  will  not,  in 
general,  be  an  exact  representation  of  the  theoretically  correct  curve; 
and  the  direction,  but  not  the  precise  location,  of  the  time  axis  of 
the  experimental  curve  is  known. 

The  generally  accepted  procedure  employed  in  analyzing  a  decay 
record  is  to  first  construct  a  theoretical  curve  that  is  a  "good" 
approximation  to  the  experimental  curve.  Then,  it  is  assumed  that 
the  damping  constants  for  the  system  are  the  damping  constants 
which  appear  in  the  set  of  parameters  that  define  the  theoretical 
curve.  Clearly,  it  is  essential  that  a  definition  of  what  constitutes 
a  good  approximation  be  given  in  terms  of  measurable  quantities.  For 
simplicity  the  goodness  of  an  approximation  is  measured  in  terms  of 
peak  displacements  on  the  theoretical  and  experimental  decay  curves. 


Statement  of  the  Problem 


With  the  preceding  discussion  as  background  the  problem  may  be 
stated  as  follows. 

With  reference  to  Figure  1,  let  XQ,  7^,  be  peak  displace¬ 

ments  on  an  experimental  decay  curve,  where  X^  Is  measured  from  an  axis 

which  Is  parallel  to,  but  not  necessarily  coincident  with  the  true 
steady  state  displacement  axis;  let  7q,  7^,  ...»  7r  be  the  corresponding 


peaks  on  a  theoretical  curve  defined  by  the  equations 

Y  -  Y  +  h  (1) 

Y  +  2vpt  +  p2tY  +  (sgn  7)A]  -  0  (2) 

Y(0)  -  Yc,  7(0)  -  0  (3) 

and,  determine  the  distance  h,  the  Initial  displacement  YQ,  and  the 
damping  constants  v  and  A  such  that  the  "square  error" 

l  n  9 

E  «  E(V,A  ,  Yq,  h)  <Xi  “  V  (4) 

l=o 


Is  a  minimum. 

To  solve  the  problem,  one  must  first  determine  how  7^  depends 

on  the  parameters  appearing  In  Equations  1,  2,  and  3.  This  Is  done 
In  the  following  section. 

Decay  Curve  Analysis 

The  extreme  values  of  the  displacement  Y(t)  which  are  obtained 
by  solving  Equation  2,  occur  at  equally  spaced  Intervals  of  time,  and 
It  Is  not  difficult  to  show  that  the  peak  displacements  are  given  by 
the  formulas 

Yt  -  Yq  (1  -  0 )  (5a) 

\  -A  (1  +  6)(—  l)1"1  -  6Yi-1  (i  +  0)  (5b) 


where 
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-  VTT 


6  •  e 


(6) 


and 

v  ■  damping  ratio  for  viscous  damping 

A  *  damping  constant  for  dry  friction  damping 

For  future  use,  it  is  noted  here  that  Equation  5  may  be  rewritten  in 
the  form 


Yt  -  Ml  +  Mai  +  Y0bt 


(5c) 


where 


o  ;  a±  -  (-l)1'1^ 


i-1 

i-i  y  6k 

kso 


b0  -  1  ;  b1  -  (-l)V 

Equations  1  and  5  may  be  combined  to  obtain 

Y  -  Y  +  h  (7a) 

o  o 


Yi  *  A(1  +  5)(-l)i_1  -  6 (?1_ 1  -  h)  +  h  (7b) 

In  view  of  Equation  6,  it  follows  that  the  relationships  expressed  by 
Equation  7  show  how  Y^  depends  on  the  damping  constants,  the  initial 

displacement  Yq>  and  h.  It  should  be  noted  that  is  independent  of 

p  (see  Equation  2)  . 

The  solution  of  Equation  7  is  of  interest.  Setting  i  *  k, 
i*k+2j,i*k+2j+l,  where  j  is  an  integer,  in  Equation  7b 
and  combining  the  resulting  equations  gives 

6  “  ”  (Yk  “  \+2j^\-l  “  7k+2j-P 
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A  “  1+V  ^*k  ”  *k+2j+l*  +  6*7k-l  ”  *k+2J^  ^ 

h  -  -  A(-l)i_1  +  (7k  +  fiYk-1)(l  +  6)"1  (10) 

and,  from  Equation  7a, 

Yo  ■  Yo  ■  h  (11) 


Equations  8  and  9  can  be  used  to  estimate  the  damping  constants  for  the 
system  characterized  by  the  experimental  decay  curve  by  substituting  X. 

r _ If  K. 


for  7, 


k* 


Y  and  h  are  Invariant  functions 
o 


Equations  8  to  11  show  that  6, A 

of  the  peak  displacements  on  the  theoretical  decay  curve.  Therefore, 
if  the  invariance  is  not  preserved  when  Y^  is  replaced  by  X^  then  it  is 


known  Chat  the  experimental  decay  curve  does  not  coincide  with  the 
theoretical  curve. 


Calculation  of  Damping  Constants  by  the  Method  of  Least  Squares 

Nothing  is  lost  and  the  subsequent  derivations  are  simplified  if 
v  is  replaced  by  6  in  Equation  4.  Then,  the  equation  becomes 


E  (6,  A,  Yq,  h)  -  (X£  -  Yp2  (12) 

The  variables  6,  4,  y  ,  h  must  satisfy  necessary  conditions  for  E  t 
be  a  minimum.  These  conditions  are , 

fi.6  -  F6(6,  A,  Yq,  h)  -  0  (13a) 

E,a-Fa(6,  A,  y0.  h)  -0  (13b) 

E,  -  F  (6,  A,  Y  ,  h)  -  0  (13c) 

xo  xo 

E,h  "  Fh(6,  V  h>  *  0  <13d> 


where  the  comma  notation  has  been  used  to  indicate  partial  derivatives 
(i.e.,  E  ,6  «8E/86,  etc.). 
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Equations  13  are  linear  in  A,  Yq  and  h,  but  nonlinear  in  6;  and 
they  are  not  amenable  to  solution  by  analytical  methods.  However,  the 
equations  can  be  solved  numerically  by  the  Newton- Raphson  Method  for 
simultaneous  equations. ^ 

To  derive  the  pertinent  equations,  and  to  outline  the  procedure, 
let  6° ,  A° ,  ,  h°  be  approximate  roots;  and,  let  d6 ,  dA,  dY^  ,  dh  be 

corrections  so  that 


6  ■  6°  +  d6 

(14a) 

A  -  A°  +  dA 

(14b) 

^  ^  +  dY0 

(14c) 

h  -  h°  +  dh 

(14d) 

Then,  Equations  13  become 

Ffi  (6°  +  d6 ,  .  .  ,  h°  +  dh)  -  0 

(15a) 

F^  (6°  +  d6 ,  .  .  ,  h°  +  dh)  »  0 

(15b) 

Fv  (6°  +  d6 ,  .  .  ,  h°  +  dh)  -  0 
*0 

(15c) 

F.  (6°  +  d6 ,  .  .  ,  h°  +  dh)  -  0 

(15d) 

Expanding  each  of  Equations  15  by  Taylor 1 s  theorem  for  a  function  of 
four  variables,  and  discarding  all  terms  containing  products  and/or 
powers  of  d6 ,  dA,  etc.,  leads  to  the  following  equations  to  be  solved 
for  the  "first  corrections,"  d6  ,  dA  ,  etc. 

*6  +  <E  66)0<l61  +  (E  A6)0dAl  +  (E  Y  6)0d^  +  (E  h6)0dhl  “  0  <16a> 

+  <E,6A)°d61  +  (E^)0^1  +  (E  ^)°d^  +  <E,hA)0dhl  "  0  <16b> 
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\  +<E,6^)°d6l+(E,AY0)0<lil+  <Ell!oYo)0d^1  +  <E,h^)0dhl  ■  0  <16‘> 

*h  +  (E,6h>°d61  +  <E,Ah>°JAl  +  +  (E,hh>°dhl  ■  0  <l6d) 

where  the  superscript  zero  indicates  that  the  quantity  is  to  be  evaluated 
at  the  "point"  (6°,  A0,  Y% ,  h°) 


Now,  the  improved  roots  are 

61  =  6°  +  dS1  (17a) 

A1  -  A0  +  dA1  (17b) 

V,1  -  Y°  +  dY,1  (17c) 

h1  -  h°  +  dh1  (17d) 


Additional  corrections  are  found  by  repeated  applications  of  Equations  13 
and  16 . 

All  possible  combinations  of  first  and  second  partial  derivatives 
of  Yi  are  needed  to  evaluate  the  differential  coefficients  appearing  in 
Equations  16.  To  calculate  the  derivatives  we  first  combine  Equations  1 
and  5c  to  obtain 

-  A(1  +  6)ai  +  Y0bt  +  h  (18) 

Then,  since  a^  and  b^  are  functions  of  6,  differentiation  of  Equation  18 
gives  (using  primes  to  indicate  differentiation  with  respect  to  6) 


?i,6  -  A[(l  +  6)a^  +  a^  +  Y0b^ 
Yi,A  "  &  +  6>ai 


(19) 

(20) 
(21) 
(22) 
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ir 

> 


7i,66  "  ^(1  +  6)®i  +  2ai^  +  Yob'i  (23) 

Yi,A6  "  +  5>ai  +  al  (24) 

Yi,Y06“bi  (25) 

All  other  second  partial  derivatives  of  are  zero. 

Evaluation  of  the  derivatives  is  simplified  by  using  recurrence 
formulas  which  follow  from  the  definitions  of  a^  and  (see  Equation  5c), 
It  is  easy  to  establish  the  following  relations 

«o  -  0  ;  a1  -  -  a±ml  +  bt  l  (26) 

bfc  -  1  5  bi  “  ~  6bi  ]L  (27) 

^  -  0  ;  a^  *  -  a^  +  b'i  ;L  (28) 

b&  «  0  ;  b*t  -  -  bi_1  -  6b't-l  (29) 


ao  0  »  ai  ai-l  +  bi-l 


b”  •  0  }  tfj  -  -  2b^j  -  (31) 

The  first  partial  derivatives  of  E  are  obtained  from  Equation  12. 
The  derivatives  are 


E,6  --S  Z1T1,6  -P6 

E,i  *-S  Zi?l,4  ‘  F4 

’  rife 

E,h--zziTl,h-*h 
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where  all  summations  are  from  i  ■  0  to  i  •  n  and  Z.  “  The 

second  partial  derivatives  appearing  in  Equations  16  are 


,66  "  2  ^  ”  ZiYi,66  +  **i,6*  ^ 

(36) 

,A6  "  2  ”  ZiYi,A6  +  Yi ,AYi  ,6^  "  E,6A 

(37) 

,Y06  2  ^  ”  Zi^i,Y06  +  ^i,Y0Yi,6^  "  E,6Y0 

(38) 

,h6  “  2  Yi,h  “  E ,6h 

(39) 

,AA  "  2  <*i,A* 

(40) 

•V  "  E  Yi,Y^7i,A  *  E,AY0 

(41) 

,hA  *  2  7i,A  =  E ,Ah 

(42) 

,YbY0  ~  2  (7i  ,Yq 

(43) 

,hY0  *  2  ri,Y0  =  E,Y0h 

(44) 

,hh  "  n  +  1 

(45) 

where  all  summations  are  from  i  ■  0  to  1  ■  n,  and  use  has  been  made  of 
the  fact  that  certain  second  partial  derivatives  of  Y^  are  identically 


zero. 


Because  of  the  special  nature  of  Equations  33,  34,  and  35,  it  is 
always  possible  to  determine  A,  Y0  and  h  as  functions  of  6 such  that 
Fa ,  F«  »  and  Fl  are  identically  zero,  as  required  by  the  necessary  con¬ 
dition^  (Equations  13)  for  E  to  be  a  minimum.  This  fact  can  be  used  to 
advantage  in  calculating  the  roots  of  Equations  13.  On  setting  F^,  F~, 

and  Fh  equal  to  zero,  combining  Equations  20,  21,  and  22  with  Equations 
33,  34,  and  35  and  rearranging  terms,  one  obtains  the  equations 


4(1  +  6)  E  a*  +  Y^  Z  a^  +  h  £  a±  -  2  J^a. 


(46) 
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A(1  +  6)  S  a^b^  +  2  b|  +  h  S  ■  L  X^b^  (^) 

A(l  +  6)Ea1+Y0  Ebi+  h(n  +  1)  -  £  ^  (48) 

where  all  summations  extend  from  1  ■  0  to  1  ■  n.  It  should  be  observed 
that  Equations  46  to  48  are  nothing  more  than  Equations  13b,  13c,  and 
13d,  but  written  in  an  alternate  form. 

Numerical  Procedure 


Equations  8  to  11  can  be  used  to  calculate  the  approximate  roots 
needed  to  start  the  numerical  solution .  However,  trial  calculations 
show  that  it  is  better  to  calculate  6  from  Equation  8  (using,  of  course, 
in  place  of  Y^) ,  and  then  calculate  A0,  and  h°by  writing  and 

solving  Equations  46  to  48  The  procedure  used  to  solve  the  problem  is 
as  follows: 

1.  Use  Equation  8  to  calculate  6°. 

2.  Calculate  A°  ,  >  and  h°  by  writing  and  solving  Equations  46  to 

48. 

3.  Calculate  Fg  using  Equation  32. 

4.  If  ■  0,  the  solution  has  been  found,  since  F^,  F^  ,  and  F^ 
are  zero  by  construction. 

5.  If  Fg  j4  0  calculate  the  differential  coefficients  in  Equations  16 
by  applying  Equations  19  to  31,  and  Equations  36  to  45. 

6.  Solve  Equations  16  for  d6*. 

7.  Replace  6° by  6°+  d6*  and  repeat  the  process. 


Damping  Ratio  and  Logarithmic  Decrement 

The  damping  ratio  and  the  logarithmic  decrement  are  measures 
of  the  amount  of  viscous  damping  in  a  system.  Let  l  be  the  logarithmic 
decrement.  It  is  shown  in  elementary  vibrations  that 


1 


(49) 


where  u  is  the  damping  ratio.  In  view  of  Equation  6  we  have 
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t  -  2Ln(6"1) 


(50) 


and,  from  Equation  49 

_  l 


(51) 


Discussion 

Sample  calculations  demonstrating  the  use  of  the  theory  are 
presented  in  Appendix  A.  In  the  example,  convergence  of  the  iterative 
procedure  was  rapid;  two  iterations  were  required.  The  decay  record 
used  in  the  example  was  taken  from  Reference  1.  Because  the  calcula¬ 
tions  are  lengthy,  only  the  first  six  (of  thirteen)  peaks  were  used  in 
the  analysis. 

The  same  decay  curve  was  analyzed  by  means  of  a  digital  computer, 
using  six  peaks  and  using  twelve  peaks.  Computer  results  are  presented 
in  Appendix  B.  A  comparison  of  the  results  Indicates  that  an  analysis 
based  on  a  partial  or  incomplete  decay  record  may  not  lead  to  accurate 
values  of  the  damping  constants. 

It  is  essential  to  note  that  an  experimental  decay  curve  may 
deviate  from  the  theoretical  curve  for  a  linear  system  with  viscous 
friction  and  dry  friction  damping  because  (a)  the  recorded  decay  curve 
is  not  an  accurate  representation  of  the  theoretical  curve;  or  (b) 
other  forms  of  damping  are  present;  or  because  of  (a)  and  (b).  The 
methods  presented  in  this  report  are  strictly  applicable  only  when 
deviation  can  be  assigned  to  cause  (a).  The  methods  may  be  applied 
when  other  forms  of  damping  are  present,  if  it  is  assumed  that  the 
calculated  values  of  6  and  A  represent  the  contribution  of  viscous 
friction  forces  and  dry  friction  forces  to  the  total  damping  force. 


CONCLUSION 

It  has  been  shown  how  the  method  of  least  squares  can  be  applied 
to  the  problem  of  analyzing  a  decay  record  to  calculate  damping  constants. 
Sample  calculations  show  that  the  method  is  not  well- suited  for  manual 
computation. 
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Appendix  A 
SAMPLE  CALCULATIONS 


In  this  appendix  the  theory  is  applied  to  the  decay  record 
shown  in  Figure  A-l.  The  decay  record  is  a  reproduction  of  the 
curve  shown  in  Figures  5-13,  page  216,  in  the  text  "Engineering 
Vibrations,"  by  L.  S.  Jacobsen  and  R.  S.  Ayre.  Because  the 
calculations  are  lengthy,  only  the  first  six  peaks  on  the  decay  curve 
were  considered.  The  procedure  given  on  page  10  was  used  in  solving 
the  problem.  The  computations  for  the  first  cycle  of  the  procedure 
are  carried  out  below.  Where  possible,  numerical  results  are 
presented  in  tabular  form. 

Calculation  of  6° 

Using  Equation  8  with  k  ■  1,  j  •  1  we  obtain 

6°  =  -  (  -  4.84  +  3. 56) (5. 60  -  4.17)'1  -  0.8951 
and ,  with  k  *  2,  j  *  1  we  have 

6°  -  -  (4.17  -  3.01)(  -  4.84  +  3.56)"1  -  0.9062 

o 

The  average  value  is  approximately  0.9,  and  6  ■  0.9  is  used  to 

start  the  solution.  The  calculation  of  a^,  b^ ,  and  their  derivatives 
is  displayed  in  Table  A~l. 
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Figure  A-l.  Decay  curve. 

(Reference:  Engineering  Vibrations,  Jacobsen  and  Ayre,  1958) 


Table  A-l.  Computation  of  a^  bi ,  and  Associated  Derivatives 


5°  -  0.9 


Column 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

Equation 

- 

- 

26 

27 

28 

29 

30 

31 

Quantity 

i 

ai 

bi 

•1 

a" 

®i 

vll 

bi 

0 

5.60 

0 

1.0 

0 

0 

0 

0 

1 

-4.84 

1.0 

-0.9 

0 

O 

H 

1 

0 

0 

2 

4.17 

-1.9 

0.81 

-1.0 

1.8 

0 

2.0 

3 

-3.56 

2.71 

-0.729 

2.8 

-2.43 

2.0 

-5.4 

4 

3.01 

-3.439 

0.6561 

-5.23 

2.916 

-7.4 

9.72 

5 

-2.42 

4.0951 

-0.5905 

8.146 

-3.2805 

17.12 

-14.58 

Equations  46  to  48  become 

40.5507  A°(l  +  6°)  -  9.0891  Y°  +  2.4661  h°-  -42.6721 

-9.0891  A°(l  +  6°)  +  3.7767  y£  +  0.2466  h°  -  19.6721 

2.4661  A°(l  +  6°)  +  0.2466  y£  +  6  h®  -  1.96000 

and  the  solution  is 

A°(l  +  6°)  *  0.2044 

^  -  5.6101 
h°  -  0.0121 
A°  -  0.1076 
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Calculation  of  First  Corrections 


The  basic  quantities  needed  to  calculate  the  first  corrections 
are  shown  in  Table  A- 2. 

Substituting  the  quantities  appearing  in  columns  5  and  10  into 
Equation  32  gives 


F l  -  —  0.1139 


and  by  the  selection  of  ^  and  h° , 


F^  -  F„°  =  F?  =  0 


Now,  Equations  16  become 

776.6387  dfi1  -  335.7554  dA1  +  41.5069  dY^ 


-  9.9599  dh  -  0.1139 


•335.7554  dfi1  +  146.3881  dA1  -  17.2692  dY*  +  4.6856  dh1  *  0 
41.5069  dfi1  -  17.2692  dA1  +  3.7767  dY<J  +  0.2466  dh1  -  0 


9.9599  dfi1  +  4.6856  dA1  +  0.2466  dY*  + 


6  dh 


And  ,  the  solution  for  d6  is 


d6  -  0.0254 


So,  the  approximate  value  of  6  to  be  used  in  the  second  cycle  of 
iteration  is 


61  *=  6°  +  d6X  -  0.9254 
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Table  A-2.  Y,  and  Associated  Partial  Derivatives 
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2.42  -2.4636  -16.2981  7.7807  -0.5905  -76.5428  19.5725  -3.2805  -0.0436 


Subsequent  Calculations 

In  the  second  cycle,  with  the  above  value  for  6*,  it  is  found  by 
writing  and  solving  Equations  46  to  48  that 

Y^  -  5.5820575 

h1  -  0.0124610 

A1  -  0.1618079 

and  F1  -  -0.002236 

6 

The  complete  solution  of  Equations  16  in  the  second  cycle  of  Iteration  is 
d62  -  0.5346  x  10'3 
dA2  -  1.1264  x  10" 3 
dY2  -  -0.5928  x  10'3 

dh2  -  0.0081  x  10' 3 

At  the  end  of  the  second  cycle  of  iteration,  the  new  approximate  value 
of  6  is 

62  -  61  +  d62  -  0.9259346 

With  the  above  value  for  6  the  solution  of  Equations  46  to  48  is 

A2  -  0.1629342 

Y2  -  5.5814647 

h2  -  0.0124691 
and  it  is  found  that 

f|  =  -0.36  x  10‘5 
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2 

Ff  is  sufficiently  small  to  justify  termination  of  the  calculation. 

o 

Table  A- 3  gives  a  summary  of  the  successive  approximations  to  the 
damping  constants. 

Table  A- 3.  Summary  of  Successive  Approximations 


Cycle 

6 

per  cent 
change 

4 

per  cent 
change 

0 

0.9 

— 

0.1076 

— 

1 

0.9254 

+2.82 

0.1618079 

+50.4 

2 

0.9259346 

40.06 

0.1629342 

+  0.7 

The  computer  program  presented  in  Appendix  B  was  used  to 
solve  the  example  problem,  and  to  analyze  the  decay  curve  shown 
in  Figure  A*1  using  the  first  twelve  peaks.  Table  A-4  shows  the 
values  for  6  as  calculated  by  hand,  and  as  calculated  with  the 
IBM  1620  digital  computer. 

Table  A-4.  Comparison  of  Calculated  Values  of  6  and  u 


No.  of  Peaks 

Manual  Solution 

Computer 

Solution 

N 

6 

8 

V 

6 

0.9259346 

0.9259393 

0.0244855 

12 

_  _  _ 

0.8952082 

0.0352147 

Consider,  now,  the  per  cent  change  in  v.  From  Table  A-4,  this 
change  is  calculated  to  be  30.5%. 
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Appendix  B 
COMPUTER  PROGRAM 


The  computer  program  was  written  in  the  Fortran  language  and 
is  shown  on  pages  21  to  24.  The  program  was  designed  to  analyse 
a  decay  record  with  a  maximum  of  30  peaks.  In  the  computer  the 

k+1  k  -  5 

iteration  procedure  is  terminated  when  d6  /  6  10  or  when 

k  ■  100  (the  maximum  number  of  iterations). 

Listings  of  typical  sets  of  input  data  are  shown  on  page  25 
and  the  output  is  shown  on  page  26  .  Printing  of  the  can  be 
suppressed  at  the  option  of  the  user  by  placing  a  zero  punch  in 
column  12  of  the  control  card.  If  the  iteration  procedure  fails 
to  converge  within  100  cycles,  the  output  consists  of  a  listing 
of  the  and  nothing  else  is  printed. 

Any  number  of  decay  records  may  be  analyzed  in  sequence,  and 
the  computer  will  halt  on  a  read  instruction  after  the  last  decay 
record  has  been  processed. 


20 


COMPUTER  PROGRAM 
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150  TYPE  160,DR,DEC,Z,N,IDENT  101 
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INPUT  DATA 


«  O  O 

w  w  w 

on  'paeo  hn^ 

*  ^U9pj  H  H  H 


*^uapx  r- 4 


fl  0% 

♦  • 

•  • 
«%  a 

►r 

A  A 

«-P  *T 

2  2 


•o  8  8 

MO  0 
CO 

0’S  ’S 

H  C8  CO 

go  a 

4J  CO  «0 
C  U  4J 
O  9  2 

OOP 


•o  o 


sumnxoQ  pa bo 


25 


OUTPUT 


CM  O 
os  CO 

<  CO  00  ,12 


CO  CM 

o  o 

I  I 

w  w 
oo  in 
Os 
r-  os 
Os  i-M  c 

co  m  |>4 


co 
o  o 

i  i 

w  w 
m  o 

CM  CM 

p  N  00  «C 

o 

I-H  CM 


CMCMCMCMCMCMCMCMCMCOCOCO 

oooooooooooo 

wwwwwwwwwwww 
ooonfONinvofiONONvo 
NflOQ  COOO^vO  o 
i-»r^o^O«nooNONOo<t<too 
|>H  CMr-vOrHi-l^t^COOOrH^r*. 
oHNnoNoooOHomn 


h  <nvo  n 
\o  oo  «h  m 
uo  <t  co 


os  so  m  cm  p  co  ^ 
«4>  o  'O  co  ON  cm  m 

N  N  H  H  C\  fs  ^ 
I  I  I  I 


CM  CM  <M  CM  CM  CM 

O  O  O  O  O  O 
U  M  M  W  M  M 

n«noooQ 

00  CO  H  01  o*  *© 

NQHgmoo 

ON  9  Q  vO  H 
CO  r-l  <o  in  CM  CO 

o>  "O  on  ^  on  n 
ino&HinoN^ 
m  ^  CO  CM  CM 
I  I  » 


26 


